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Results are presented for spin-wave dispersions in geometrically frustrated stacked triangular 
antiferromagnets with a thin film or semi-infinite geometry having either zero, easy-plane, or easy- 
axis anisotropy. Surface effects on the equilibrium spin configurations and excitation spectrum are 
investigated for the case of antiferromagnetically coupled films, serving to extend previous results on 
ferromagnetically coupled layers [E. Meloche et ai, Phys. Rev. B 74, 094424 (2006)]. An operator 
equation of motion formalism is applied to systems which are quasi-one and quasi-two dimensional 
in character. In contrast to the case of ferromagnetically coupled films the new results show surface 
modes that are well separated in frequency from bulk excitations. Magnetic excitations in thin films 
with an even or an odd number of layers show qualitatively different behavior. These results are 
relevant for a wide variety of stacked triangular antiferromagnetics materials. 

PACS numbers: 75.30.Ds, 75.70.-i 



I. INTRODUCTION 



The significant amount of research effort devoted to the discovery and understanding of new aspects of geometrically 
frustrated magnetic systems in recent decades has been almost exclusively focused on the study of bulk propertiesii 
Among such systems, a large number of materials have been identified which are realizations of the prototype frustrated 
triangular antiferromagnet (AF). In addition to the AF near-neighbor in- plane coupling, the vast majority of these 
systems are also characterized by AF coupling between stacked triangular layers, such as in most of the ABX3 
compounds.^ Theoretical and experimental studies of magnetic excitations in these model materials have revealed 
important information on the nature of the fundamental interactions and exposed new types of spin-wave modes in 
the bulk systemsi^i^'^'^ Interest in surface effects on frustrated magnetic systems has been enhanced recently with 
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the fabrication of new materials classes. Studies of hexagonal magnetoelectric and multiferroic rare-earth magnetites 
(e.g., RMnOa) thin films are being fueled by potential technological applications in the field of spintronics. - Ultrathin 
films of Mn on fee (111) Fc-Ni substrates also offer the possibility to study AF triangular layersj^ The understanding 
of thin-film, surface and interface spin dynamics of such systems is of interest for the engineering of high data transfer 
rate devices<^ 

In the present work, results arc shown of solving the operator equation of motion to examine surface spin excitations 
in a model of AF coupled frustrated triangular layers which have either easy-plane or easy-axis anisotropy. This 
study represents an extension of our previous work where the technique was developed and applied to the case 
ferromagnetically (F) coupled triangular layers having easy-plane anisotropy [E. Meloche et ai, Phys. Rev. B 74, 
094424 (2006)]. (A summary of general theoretical and experimental techniques that have been employed to study 
thin-film and surface spin waves in un-frustrated systems is also provided in Ref. [10].) The formulation of a linear 
spin wave theory is based on a model Hamiltonian which includes isotropic and anisotropic exchange as well as single- 
ion anisotropy and allows for layer-dependent coefhcients of these interactions. Depending on the specific system of 
interest (e.g., semi-infinite and thin films with or without adjacent layers of a different material) surface parameters 
may have values quite different from the bulk. A significant difference found for the present case of AF-coupled layers 
is that localized modes can be well separated in frequency even with bulk parameter values used at the surface layers, 
in contrast with F-couplcd films. Also, as might be expected, equilibrium spin configurations and excitations are 
dependent on whether the films are composed of an even or odd number of layers (again in contrast with F-coupled 
films). Due to physical geometrical considerations, nearly all material realizations of stacked triangular layers can be 
classified as having either quasi-one or quasi-two dimensional magnetic character. Illustrative numerical results of the 
spin wave calculations for these cases are discussed in the following sections with the intent to reveal key aspects of 
thin-film and surface effects on spin excitations in these systems. 

The systems considered here may be represented using the following spin Hamiltonian 

n= J^^i {^^ ■ ^0 + ^^ts^) + • + '^'stst) + A(5n' (1) 

<i,j> <i,i'> i 

where Jij > represents the intralayer nearest-neighbor AF exchange coupling and J'- j/ > represents the interlayer 
AF exchange coupling. The effects of anisotropic exchange coupling are included through the parameters a and cr', 
and Di represents the strength of the single-ion anisotropy at a magnetic site labelled i. The Hamiltonian in Eq. ([T]) 
can used to describe systems characterized with easy-plane {Di > 0) and easy-axis {Di < 0) single-ion anisotropy. 



Easy-plane and easy axis systems will have different equilibrium spin configurations and spin dynamics and therefore 
the cases will be treated separately. In Sec. HH bulk and surface spin-wave modes are calculated for the case of 
easy-plane anisotropy. In Sec. IIIII results are deduced for the spin-wave energies and spectral intensities for the case 
of easy-axis anisotropy. A discussion of these results and their relevance to specific compounds is presented in Sec. 
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II. EASY-PLANE ANISOTROPY 
A. Easy-plane bulk spin-wave excitations 

Extending the general method described in Ref. [11] to calculate the bulk spin- wave dispersion relation, the system 
is divided into three pairs of sublattices on two adjacent layers, which are labeled Ai, Bi, Ci; A2, B2, and C2, and the 
specific sites on each sublattice are labeled with indices Zi, mi, rii, I2, TO2, and n2, respectively. Within any particular 
layer the 120° structure is stabilized as a result of the pairwise antiferromagnetic coupling between nearest- neighbor 
spins. Also, in the classical ground state configuration the spins on adjacent layers (along the crystallographic c 
axis) order antiparallel to one another. We first transform to a local coordinate system such that the z-axis for each 
sublattice is in the direction of the average spin alignment. The transformations to the local coordinate systems for 
the various sites may be written as 

^ [{-ir+\sf cos 0,- Sf sin 0i),{~ir+\sf sine, + Sf cos 9,),-Sf'] (2) 

with i ~ l^^mf^fif^ and = 0, — 27r/3, — — 27r/3 with /i = 1,2. The linearized equations of motion for 
the local spin raising and lowering operators at each lattice site are formed and then transformed to a wavevector 
representation. Assuming the usual time dependence exp(—iEt), the resulting set of 12 equations may be expressed 
as Mb = 0, where 



and b — 



bi 

b2 

The elements of the 6x6 block-circulant matrices A and t are written as 




(3) 
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^ A B B*^ 



A = 



B* A B 
B B* A 



^ T{k,) ^ 

T(fc,) 
\ T(fc^) J 



The matrix elements of A and r depend on the energy and system parameters and are defined as 



E + n 






% -A 






, B = 




E-Q. J 




^-7 -Hj 



and T{k^)=SJ'{k^)/2 



-a' 2 + a' 
-2 - a' a' 



where 



Q. = -SJ{0)- SD' - SJ'{0) 

a = SD' 

P = -5(1/2 + a)J(k||)/2 

7 = 5(3/2 + a)J(k||)/2 

D' = [1-{2S)-^]D. 



The exchange integrals J(k||) and J'{kz) are defined as 



(4) 



(5) 



(6) 



J(k|| ) = J(2 cos(fcy \/3a/2) exp(— ifcxa/2) + exp(ifcxa)) 
J'{k^) = 2J'cos{k^c) 



(7) 



with a and c denoting the lattice constants. The elements of the column vector b are written as 



= [sX, (k), si^ (-k), s+ (k), 5^., (-k), s+ (k), (~k)] 



(8) 



for /i = 1,2. The terms S'^^(k) correspond to the Fourier amplitudes of the spin operators S^, for /x = 1,2, along 
with similar definitions for S'^^(k) and S'^^(k). The bulk spin- wave modes correspond to the solutions of detM = 0. 
Analytical results may be obtained by first block diagonalizing M using the transformation W^MW — D where 



W= — 

V2 



u u 
-u u 



(9) 
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U = 



1 

7! 



with X = cxp(i27r/3) 



(10) 



I2 I2 I2 

\ I2 X*t2 Xt2 J 

The nonzero elements of D consist of six 2x2 matrices along the main diagonal and the bulk spin-wave energies are 
obtained from the solutions of the determinantal conditions det(A(0) ±t) — for (p — 0, ±2tt/3 where we have defined 
A{ip) = A + B cxp(i(/)) + B* exp(— Bulk spin-wave energies may be expressed as Wk±(0) with = 0, ±27r/3, where 



tul±{(j)) = {fl + a± SJ'ih) + S'J$(k||, (/))) {n-aTS{l + a')J'{h) - 2SJ{1 + CT)*(k|| , 0)) 



(11) 



and 



$(k||,0) = 2coa{V3kya/2) cos{-k^a/2 + (/)) + cos{k.^a + (f)) . (12) 

B. Easy-plane surface spin- waves 

As in our previous work,^^ calculations are first made for a semi-infinite system with a single (001) surface and 
for a thin film composed of N magnetic layers which are labeled using the layer index n (= 1, N). The linearized 
equations of motion for the local spin raising and lowering operators at all the sites are formed and then transformed 
to a representation involving a two-dimensional wavevector k|| which runs parallel to the surfaces and a layer index 
n. The spins located on the layers with an odd layer index n are taken to belong to sublattices Ai, Bi, Ci whereas 
those located on layers with an even layer index n belong to sublattices A2, i?2, ^2- For the l^^ site on sublattice Ai 
the wave-like solution for the spin operator is written as 



S^,=SX^^{}^\\)eMi{^WP-m] (13) 

where the position vector p = {x, y) and the amplitudes S*^^ nC'll) depend on the z coordinate through the layer index 
n. Similar expressions are defined for sites on sublattices Bi,Ci, A2, B2 and C2. For a semi- infinite system the set 
of finite difference equations connecting the Fourier amplitudes may be expressed in supermatrix form as Mb = 0, 
where M is an cx) x 00 block-tridiagonal matrix defined as 
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M = 



Ai t(0)/2 
t(0)/2 a t(0)/2 
t(0)/2 a r(0)/2 



V 



(14) 



where each element M represents a 6 x 6 matrix and A is defined as in (|4|)-(l7|). The elements of the matrix Ai are 
defined as in A but with the substitution J ^ Ji, D ^ Di and J' J' /2. The differences between Ai and A occur 
because of different exchange and anisotropy parameters for spins in the surface layer as well as the reduced number 
interlayer neighbors for surface spins. The infinite column vector b is written as b = [bi,b2, ...]"^ where b„ is defined 
in Eq. ([H]) and n — 1, oo. The equations are first partially decoupled by applying the transformation U^^Mi.jU to 
all of the elements of the supermatrix. The transformed equations may be written into three sets of finite difference 
equations as 



Ai(0)Xi,^ + ^X2,^ = {n^l) (15) 
^X2.,^+A(0)X2.+i,^ + ^X2,+2.0 = (n-2z.+ l,i.>l) (16) 

^X2.-i.0+A(,/))X2,,^ + ^X2.+i,^ = {n = 2iy,iy>l) (17) 

for = 0, ±27r/3. The matrix Ai{4>) is defined as Ai(0) = Ai + Bi exp(i0) + exp{—i(j)), whereas A(0) is defined as 
in the bulk case. The column vectors in Eq. p5|) - p?|) are written as X„_0 — [ AT +((/>), AT^ (</>)] with elements 



X±(0) = ^(^i„ + ^|,„ expi-z^) + 5±„ exp(z0)). (18) 

where we omit the wavevector label k|| . Column vectors X„.0 with an odd (even) layer index n involve linear 
combinations of the operator amplitudes on sublattice 1 (2) only. Equations and (|16p are used to eliminate the 
^n,4> with n odd from the set of equations represented by (|17p . The resulting set of relations connecting X„^0 with 
n even can be written in the matrix form 
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/ 



Pi,0 Q0 
Q0 P0 Q0 
Q(f> P4, 



V 



X 



4,0 



(19) 



V ■ J 



where 



P1.0 = A(0)--r(O)(Ar^(0)+A-i(0))T(O) 
P0 = A(0)-ir(O)A-i(0)r(O) 

= -ir(O)A-i(0)r(O) (20) 

The three sets of equations given in Eq (fT9|) (for (/) = 0, ±27r/3) couple the amphtudes on sublattice 2 only. The surface 
spin-wave frequencies can be obtained numerically following an analogous approach employed for ferromagnetically 
coupled layers.— 

For a thin film composed of N layers with (001) surfaces the system of finite difference equations can be expressed 
in a supermatrix form as Mb = as for the semi-infite case (see Eq. [H]) except M is now a 6iV x 6iV block-tridiagonal 
matrix and the column vector is defined as b = [bi, ...jh^]"^. The first and the last elements on the main diagonal 
of the supermatrix M are written as Ai and Ajy, respectively. The elements of the matrix An describe the effects 
of the additional surface of the film and are defined as in Ai for the semi-infinite case except with the appropriate 
exchange and anisotropy parameters at the surface n = N . The coupled systems of equations in the case of a thin 
film can be decoupled using the method outlined for a semi-infinite system. However, cases involving an even or an 
odd number of layers N must be dealt with seperately. For a film containing an odd number of layers both surfaces 
belong to the same sublattice whereas for a film an even number of layers the spins in the surface layers belong to 
difference sublattices. The spin-wave excitation spectrum in a stacked triangular antifcrromagnctic film are obtained 
numerically by solving the determinantal condition det M = 0. 



C. Numerical Results 

In Fig. [T] illustrative results are shown for the spin-wave frequency versus in-plane wavevector kxO, (with ky ~ 0) 
for a representative S = 1 quasi-lD planar triangular antiferromagnet where we set J = 1.0 GHz, J' — 100.0 GHz 
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and (T = a' = 0.0. The dispersion relations in the top and bottom figures correspond to D' = 1.0 GHz and D' = 0, 
respectively. The surface exchange and anisotropy parameters are taken to be uniform throughout the system. The 
shaded areas represent all of the bulk solutions a;k±(0) (see Eq. (fTTj) ) with < kzC < 7r/2 and the three sections 
(from left to right) correspond to the bulk solutions with = 0, 27r/3, — 27r/3. The Af-point denotes the zone-edge 
wavevector {kx,ky) — (27r/3a,0). For these parameter values the top of the bulk regions extends to approximately 
208 GHz and only the lower part of the effective bulk continuum is shown. The solid line represents a surface spin- 
wave obtained from the solution of Eq. (jl9[) for a semi-infinite system. This new mode is characterized with an 
amplitude that decays into the film away from the surface. In the top figure (zero anisotropy) results are also shown 
for two lowest energy modes for thin films composed layers 8 (dotted lines) and 16 layers (dashed lines). In the 
bottom figure the dotted, dashed and solid lines represents two lowest energy excitations for a film composed N —9, 
17 and 100 layers, respectively. The dominant effect of the anisotropy is to open up a gap in the spectrum at the 
F-point for the surface and bulk excitation. The gap energies for the bulk modes are obtained from Eq. (jlip and are 
cjk-(±27r/3) = S[D'{8J' + 9J)]^^^. For most values of the wavevector the surface mode obtained for a semi-infinite 
system is well separated from the bulk region. This contrasts with the behavior obtained for a quasi-lD easy-plane 
triangular antiferromagnet with ferromagnetic interlayer exchange coupling where the splitting between the surface 
branch and the bulk region is negligible for a system characterized with uniform exchange and anisotropy parameters. 

The spin-wave dispersion relation in these frustrated films show interesting characteristics. The dispersion relation 
in films composed of odd number of layers are qualitatively similar to those with ferromagnetic interlayer exchange. 
This can be explained by the fact that for an odd number of layers both surfaces belong to the same sublattice and 
each individual chain of spins will possess a net moment. For thin films with an even number of layers N the system 
is composed of antiferromagnetically coupled chains that possess no net moment because each chain of spins has an 
equal number of sites from sublattices 1 and 2. As the number of layers in the film increases this effect becomes 
less important and the two lowest branches eventually become degenerate and equal to the results obtained for the 
semi-infinite system. The two modes correspond to localized excitations at each surface. 

In Fig. [5] results are shown for the surface and bulk spin-wave energies versus in-plane wavevector k^a (with 
ky = 0) for a semi-infinite quasi-2D system with J = 100.0 GHz, J' = 1.0 GHz, D' = 1.0 GHz (with Di = D) and 
a = a' = 0.0. The bulk spin waves form a very narrow continuum (shaded areas) in quasi-2D systems because of the 
weak dependence on the third wavevector component kz ■ Also shown are the substantial effects of a modified surface 
exchange coupling Ji on the surface spin waves. Splitting of the surface branch away from the bulk region can also 
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FIG. 1: Spin-wave energy versus in-plane wavevector kxa for a quasi-lD system with J = 1.0 GHz, J' = 100.0 GHz and 
a = a' = 0.0 and in the top figure D' = 1.0 GHz whereas in the bottom D' = 0.0. The M-point refers to the wavevector 
k|| = (27r/3, 0). Shaded regions correspond to the bulk excitations. In the top figure, solid, dotted and dashed curves correspond 
to the semi-infinite, 8-layer and 16 layers systems, respectively. In the bottom figure, dotted, dashed and solid lines correspond 
to films with 9,17 and 100 layers, respectively. 

be obtained by assuming a modified value for tlie surface anisotropy Di compared to the bulk value. The energy gap 
at the zone center F for the bulk modes a;k-(±27r/3) is due to the anisotropy and vanishes when D' = 0. 
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FIG. 2: Surface and bulk spin-wave energy versus in-plane wavevector k^a for a semi-infinite quasi-2D system with J = 100.0 
GHz, J' = 1.0 GHz, and D' — 1.0 GHz (with Di = D) and a = a' = 0.0. The solid line appearing just below the bulk region 
(shaded area) corresponds to a case with uniform interlayer coupling Ji = J, whereas the dotted and dashed lines are obtained 
with Ji = 0.5J and Ji = 0.75J, respectively. 

III. EASY-AXIS ANISOTROPY 
A. Easy-axis bulk spin-wave excitations 

Consider now the triangular antiferromagnet characterized with a singlc-ion anisotropy whose easy-axis is along 
the crystallographic c axis. It is instructive to first carry out the calculations for the bulk spin-wave dispersion. We 
follow a method similar to the one described previously for the case of easy-plane anisotropy. For the bulk spin- wave 
calculation, the system can again be divided into the six sublattices as depicted schematically in Fig. [3^. It is worth 
noting that the ground state spin configuration can be rotated about the c axis by any arbitrary angle and here we 
assume that all of the spins have their equilibrium positions in the ac plane. Without loss of generality we set the 
equilibrium direction for spins on sublattice Ai along the c axis. The equilibrium configurations for spins on sublattices 
Bi,Ci, A2, B2,C2 can be obtained using a local field method.i^ii^ The main effect of the easy-axis anisotropy on the 
equilibrium configuration is that the spins now lie in the a plane containing the c axis and a deformation of the 120° 
structure giving 6 > 120° as shown in Fig. [3)3. 

The Hamiltonian in Eq. ^ is transformed to a local coordinate system such that the direction of the local z-axis is 
in the direction of the equilibrium spin alignment. The transformations for the sites on the six sublattices are written 
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FIG. 3; (a) The deformed 120° structure for a STA with easy-axis anisotropy. The angle 6 is greater than 27r/3 when D < 
(when a = g' = 0.0). AU of the spin he in the ac plane and the direction of the spins on sublattices A2, B2, C2 are opposite 
to those on sublattice Ai, Bi, C'l, respectively. 

as 

{Sf,Sf,Sn -> [{-ir+\sf COS 0,- Sf sine,), Sf ,{-ir+\sf sin 0, + Sf cose,)] (21) 

with i = /^t, m^, and 0i^ — 0, 0m^ = 0, = —d with ^ — 1,2. The stability condition requires that the coefficients 
of the local transverse spin components Sf and Sf vanish for all sites and leads to the bulk canting angle which is 
given by 



^° 2[1 + a/2 - a'J'/'SJ + D'/3J] ' ^ ' 

Eq. (|22p is a generalization of the results obtained in Refs. [13,14] which treat the cases of either easy-axis anisotropy 
or exchange anisotropy, but not both. In the absence of anisotropy {a = a' = D' = 0) the 120° structure is recovered, 
as expected. The linearized equations of motion for the various sites are calculated and then transformed to a 
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wavevector representation. Twelve equations are required to obtained a closed set. The system of coupled equations 
may be expressed as Mb = where the 12 x 12 matrix M has the same form as in Eq. ([3]) and the column vector b 
is formally written as in Eq. ([8]) . In this case, the elements of the supermatrix M are now given by 



B* C D 
B D* C 



^ T{k,) ^ 

A(fc,) 
\ A(fc,) j 



where 



A 



C 



T{k,) = Sj'ik 




B = -5J(k||)/2 



D 



-5J(k||)/2 



A(fc,) = SJ'{k,)/2 





:;) 








-:;) 




a'sin^ 



-{2 + a' sin^ 



(2 + cr' sin^ 
-a' sin^ ( 



with matrix elements defined as 

n 



s = 



2S{1 + a) cos6lj(0) + 2SD' - 5(1 + ct') J'(0) 

^((l + ct)(cos6' + cos^ 0) - sin^ 0)J{O) - 5(1 - 3 cos^ e)D' - 5(1 + ct' cos^ 0)J'{O) 
cos 6* ± 1 

cos^ 6* - (1 + ct) sin^ 9 ± 1 
D' = [1-{2S)-^]D. 



-1 



(23) 



(24) 



(25) 



The exchange integrals J(k||) and J'{kz) are defined as in Eq. ([7]). The matrix A is not a block-circulant matrix 
due to the effects of the anisotropy and analytical solutions for the bulk spin-waves cannot be obtained using the 
diagonalization procedure employed in Eq. (Pl)- (|10p . Instead, the bulk spin- wave spectrum is obtained numerically 
by solving the determinantal condition detM = 0. 
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FIG. 4: (a) Schematic view of the spin configuration for a stacked triangular antiferromagnet with easy-axis anisotropy and a 
(001) surface as viewed along the chains, (b) Local coordinate system for spins in layer n, where n is an odd layer index. 



B. Spin-waves in thin films 

The spin-wave excitations in thin films composed of magnetic layers with a (001) surfaces may be investigated 
by extending the method used for bulk spin-waves. The calculations for a thin film with easy-axis anisotropy are 
algebraically more complicated than those for a film with easy-plane anisotropy because the system can no longer be 
described in terms of a six sublattice model. However, within any particular layer of the film the system is still divided 
into three interpenetrating sublattices. Without loss of generality, the equilibrium direction for the surface layer spins 
n = 1 on sublattice A is set along the c axis and the equilibrium direction for all of the other spins are obtained 
by miminizing the classical energy of the system. Sites on sublattice A with an odd layer index n also have their 
equilibrium directions along the c axis, whereas sublattice A sites with an even layer index are antiparallel because 
of the antiferromagnetic interlayer exchange coupling. The canting angles for sites on sublattices B and C in layer n 
are written as Ob.u = dc,n = dn- A schematic diagram of the system is illustrated in Fig. 0^ along with the local set 
of axes for each sublattice. The 3A^ coordinate transformations may be written as 
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for i = lmm„,nn with n — 1, ...,-/V, and we define = 0, 0m„ = ^n, ^n„ = ~^n- The Hamiltonian in Eq. ([1]) is 
again written with respect to the local coordinate systems and the stability conditions for the spins on sublattices B 
(or C) are now given by 



3J„[(l + cr) sin6l„ + (l + (7/2) sin 26'n] + D'^ sin 26'„ 
+J' [(sin0„_i cos0„ — (1 + (t') cos0„_i sin0„)(l — &\.n) 
(sin^n+i co^Bri - (1 + cr')cos6'„+i sin6'„)(l - (5„,jv)] = 

for n=l,...,7V. 



(27) 



The surface effects are taken into account through the Kronecker 6 in Eq. (P7|) and modified values of the exchange 
and single-ion anisotropy at the surfaces of the film are allowed. The set of N nonlinear equations can be solved using 
a numerical self-consistent approach. ^^'^^ However, there are a few limiting cases that can be solved analytically. If 
we neglect the anisotropy terms (a = a' = D'^ = 0) we recover the 120° structure in any particular layer, as expected. 
In the limit of J' = the layers are independent from one another and the canting angle within any layer will be 
given by Eq. (|22p . When the exchange and anisotropy surface parameters are set equal to their bulk values it is easily 
verified that the canting angles for sites on sublattices B and C are layer independent and again given by Eq. (|22p . 
For all other numerical solution is required. 

The spin- wave energies in thin films are obtained by forming 6-/V linearized coupled equations of motion which may 
be expressed as Mb = where M is a block tridiagonal matrix defined as 



M = 



Ai ri,2 

T\,1 A2 'r2,3 







V 



TiV-2,JV-l Aat-i rN-l,N 
• • • TN-1,N Aat 



(28) 



and b = [bi, ...,bjv]"^. The elements of the matrices A„ and x„_„/ (defined in Appendix A) are similar to those in 
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FIG. 5; Canting angle (measured in degrees) for sites on sublattices B or C versus layer index for a A'' = 25 layers film with 
easy-axis anisotropy. All the curves correspond to D'l = D' — D'j^ — —J and a = a' = 0. The solid line corresponds to the 
case Ji — J — Jn and the canting angle (Eq. (I22|l 'l is independent of the coupling strength J'. The dashed line and the circles 
show results for the quasi-lD cases of J' = lOOJ and J' — lOOOJ when the surface exchange couplings differ slightly from the 
bulk values (Ji = Jm = 0.9 J). 

the bulk case in Eqs. (j23p - (|25p where the main differences are due to the layer dependent canting angles for spins on 
sublattices B and C. Spin-wave excitations are obtained numerically from the solutions of det M = 0. 

C. Numerical Results 

In Fig. [5] results are shown for the canting angle for spins occupying sites on sublattice i? or C in a 25-layer film 
with easy-axis anisotropy. In all cases the spins on sublattice A form antiferromagnetic chains with ordering along 
the c axis. The canting angle is layer dependent only if the surface exchange or anisotropy parameters are allowed 
to differ from their corresponding bulk values. In quasi- ID systems the strong interlayer coupling tends to restore 
the antiferromagnetic order in the spin chains on sublattices B and C. For quasi-2D systems the perturbation of the 
canting angle in cases with modified surface parameters is localized to a small number of layers (typically ~ 2) away 
from the surface. 

An example of the spin-wave dispersion relation for a quasi-lD system with easy-axis anisotropy is illustrated in 
Fig. [6l The shaded area corresponds to the low-energy bulk spin-waves. The solid lines represent the six lowest 



16 




FIG. 6: Spin-wave energy versus in-plane wavevector feo (with ky = 0) for a S = 1 quasi-lD system with easy-axis anisotropy. 
The parameters are J = 100.0 GHz, J' = 1.0 GHz, D' = —1.0 GHz, a = a' = 0.0. The exchange and anisotropy parameters 
at the surfaces of the film are set equal to their corresponding bulk values. The shaded area corresponds to the lower edge of 
the bulk continuum and the solid linos correspond the six lowest energy modes obtained for a thin film composed of Af = 20 
layers. The six lowest energy branches for the film are labelled u}i,...,u)e. 

energy modes for a uniform film composed of A'' = 20 layers. The parameter values for the representative quasi-lD 
system are J = 1.0 GHz, J' = 100.0 GHz, D' = —1.0 GHz and a = a' = 0.0. All of the spin-waves modes obtained 
from the determinental condition det M = are degenerate in magnitude and only the positive solutions arc shown. 

The operator equation of motion method used in the previous sections may be easily extended to calculate Green 
functions which can be used to discuss the spectral intensities of the new spin-waves modes predicted in this work. 
In order to further characterize the low-energy spin-wave modes in thin films with easy-axis anisotropy the mean- 
squared amplitude of the spin precession, defined for sublattice i (= A,B,C) as (5*(k||,n) = {{Sf^^)^ + ('S'f„)^)k|| , 
are considered. For simplicity, we outline the calculation of (5'^(k||,n). This quantity can be written in terms of the 
transverse spin-correlation functions {S^ ^{t)S^ ^{t'))k^^ and (5'J„(t)5'^„(i'))k|| evaluated at equal times {t = t'). 
The spectral representation ^n(k|| , E) of the transverse spin-correlation function {S^ ni^)^A n('^'))^i\ defined as 

/oo 
Uh,E)exp[-iE{t-t')]dE, (29) 
-OO 

is obtained from the imaginary part of the Green frmction Gj^^(kj|, E) =<C S^rii^w)'-' '^a n(^ll) using the fluctua- 
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FIG. 7: Mean-squared amplitudes Q^(k||,n) = Q'~'(k|[,n) plotted against the layer index n for the some of the low-energy 
k|| = (27r/3a, 0) modes defined in Fig. |6]for a quasi-lD film composed of A'^ = 20 layers. The amplitudes (3^(k||,n) of the 
spin- wave modes UJ2 and los (not shown) are qualitatively similar to those of uji . All of the results are symmetric with respect 
to the center of the film because of the symmetric choice of anisotropy and exchange parameters at the surfaces. The lines are 
guides to the eye. The closed symbols define the amplitudes on even numbered layers whereas the open symbols are used for 
the odd numbered layers. 

tion dissipation theorem.— The spectral function of (S"^ „(i)'S'^ „(i'))k|| may also be deduced from G;^„(k|| , E). Some 
examples of the application of the Green function cquation-of-motion method to magnetic thin films with a different 
ordering of the spins are given in Refs. [18-20]. In the low-temperature limit {T Tc) the solution to the Green 
function G^„(kj|,i?) may be written as 



G;^_„(k||,£;) = -^-^[adjM]6„-5.6n-5 for n=l,...,7V (30) 

where [adj M]„.„ denotes the n**^ diagonal element of the adjoint of the matrix M defined in Eq. ([28|l . Other diagonal 
elements of the adjoint of M are used to obtain the Green functions G^„(k||,i?) and G5^„(k|| , i?). The spin-wave 
energies obtained from det M = will correspond to the poles of the Green functions. The mean-squared amplitude 
(5"^(k|| , n) (or integrated intensities) of a particular spin- wave mode is estimated by evaluating the area under the peak 
of the imaginary parts of G^„(k||,£' + ie), where the real and positive quantity e is introduced phenomenologically 
to model an intrinsic damping or reciprocal lifetime. 

Fig. [7] shows the mean squared amplitudes as a function of the layer index corresponding to some of the k|| = 
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(27r/3a, 0) spin-wave modes shown in Fig. [6l The mean-squared amplitudes for sublattices B and C within any 
particular layer are equal. Numerical results show that the spin- wave modes loi, ...,a;4 contribute mostly to the spin 
fluctuations on sublattices B and C. These modes are localized near the surfaces of the film and are characterized 
with amplitudes which decay with distance from the surfaces. The amplitudes Q^(k^^,n) (not shown) of the spin- 
wave modes wi, ...,104 are negligible. The differences in amplitudes between odd and even layers is a result of the AF 
ordering of the spin chains. The spin- wave modes W5 and ujq contribute much less to the fluctuations on sublattices 
B or C compared with modes wi, ...,0^4. These higher energy modes are characterized with amplitudes that vary in 
a wavelike fashion across the thickness of the film on sublattices B and C. However, numerical results show that 
the spin-wave modes 0^5 and wg are localized on sublattice A sites that are near the surfaces the film. Higher energy 
modes are characterized with amplitudes that vary in a wavelike fashion on all three sublattices. 

Fig. [H] shows illustrative results for a 5 = 1 quasi-2D film composed of = 20 layers where we set D' = —10 
GHz, J' — 1.0 GHz, J = 100.0 GHz <j = a' = 0.0. Exchange and anisotropy parameters at the surfaces are set equal 
to their corresponding bulk values. For each value of the in-plane wavevector k|| the solution of det M = yields a 
total of 6iV spin-wave modes. All of the solutions are degenerate in magnitude and the dispersion relation (showing 
positive solutions only) is illustrated in Fig. [51A_. The modes are split into 3 distinct groups each containing a total of 
20 branches. Qualitatively similar results (not shown) are obtained for a quasi-2D system with easy-plane anisotropy. 
The main effect of the easy-axis anisotropy is to remove the degeneracies between some of the spin modes. This is 
a consequence of the reduced symmetry in systems characterized with easy-axis anisotropy compared with systems 
with easy-plane anisotropy. Discussions of the symmetry properties of the bulk spin-wave modes can be found in 
Ref. [13]. The lowest energy branch in each group is labelled wi, W2, and, ^3. In quasi-2D systems these modes 
are two-fold degenerate and each mode represents a localized excitation at the surfaces of the film. The 18 other 
spin-waves branches within each group appear in the effective bulk continuum (not shown) obtained for the infinite 
system. The inset in Fig. shows the splitting in frequency (at small wavevectors) of the surface branch labelled 
LU3 from the rest of the modes in the group. In Fig. [8}d results are shown for the spectral representation obtained 
from — Im{G;^„(k|| , i? + ie)} for the 20 low-lying k|| = (27r/3a,0) spin-wave excitations around uji. The solid and 
dashed lines show results for the spectral representation in layers n = 1 and n = 3, respectively. The intensity of 
the spin-wave mode uji is largest in the surface layers (n = 1 or n = A^) and decays extremely rapidly with distance 
away from the surfaces. In Fig. [8j: we show the amplitude (5'^(k|| , n) (integrated intensity) of the surface modes as a 
function of the distance away from the surface. The results are symmetric about the center of the film and only the 
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Frequency (GHz) 




FIG. 8: (a): Spin- wave dispersion relation for sl S — 1 easy-axis quasi-2D film composed of 20 layers with uniform parameters 
where we set _D' = — 10 GHz, J' = 1.0 GHz, J — 100.0 GHz, and a = a' = 0.0. The spin- wave modes are split into 3 groups 
each containing 20 branches. The labels uji, UJ2, ^3 refer to the lowest energy branches in each group. The inset focuses on 
the small wavevector region and shows the splitting of the uj-i, branch from the other modes in the group, (b): Sublattice A 
spectral representation (— Im{G'^_,j(k|| , i5 + ie)} with ky = (27r/3a,0) and e = 0.01 GHz) versus frequency in layers n — 1 (solid 
line) and n = 3 (dashed line), (c): The amplitude Q'*(k||,n) (integrated intensity) of the zone-edge surface modes uji, lu2, UJ3 
as a function of the layer index. Results are symmetric about the center of the film and only results from the first 4 layers are 
shown. Lines are guides to the eye and have no physical meaning, (d): The mean-squared amplitude (3^(k||,n) for the bulk 
excitation ljb defined in Fig. |8]d versus layer index. Open and close symbols are used for odd and even layers, respectively. 
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FIG. 9; (a): As in Fig. |8] except with modified exchange parameters at the surface of the film. Results correspond to 
Ji ~ Jn ~ 0.9 J. The labels uji, ui2, refer to the 3 surface modes, (b): Sublattice A mean-square amplitude evaluated at 
k|| — (27r/3a, 0) for the surface modes as a function of layer index showing the first 4 layers only. The amplitudes are ~ 3 
orders of magnitude smaller than those in Fig. [8};. 

first 4 layers are shown. Qualitatively similar behavior is obtained for the amplitudes Q^(k||,n) (= Q'-^(k|| , n)). The 
surface modes have negligible spectral intensities for any layer 4 < n < 16. Fig[8]i shows results for the amplitude 
as a function of the layer index for the bulk mode labelled lob of Fig. [8}d. The modes appearing within the effective 
bulk continuum have amplitudes that vary in a wavelike fashion across the thickness of the film. The intensities of 
the bulk excitations at the surfaces is strongly dependent on wavevector and the assumed anisotropy and exchange 
parameters at the surfaces of the film. 

Some effects of modified values of the surface exchange parameters on the spin-wave dispersion relation and the 
amplitude of the modes are illustrated in Fig. [9] for a quasi-2D system. The results are obtained using the same 
parameters as in Fig. [5^ except the exchange parameters at the surface are set as Ji — Jn — 0.9J. The localized 
surface modes labelled lui, lu2, and, W3 are shifted to lower frequencies when the surface exchange parameters are 
less than the bulk values. Fig. shows that the modes uji, uj2, and uj^, are localized near the surfaces of the film 
and are characterized with amplitudes that are orders of magnitude smaller than the case with uniform parameters 
throughout the thickness of the film. For a symmetric film with Ji = Jn and Di = Dn the surface modes are two-fold 
degenerate. However, the degeneracy is lifted when the surfaces are asymmetric. 
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IV. CONCLUSIONS 

The results presented above outline the application of an operator formalism to calculate of surface and thin-film 
effects on the spin dynamics in AF coupled geometrically frustrated triangular layers. This serves to extend our 
previously published work on F coupled films to the physically more relevant case applicable to a large variety of 
materials. Through the use of illustrative numerical calculations differences in linear spin excitations for bulk, semi- 
infinite and thin-film systems having either easy-plane or easy-axis anisotropy and quasi-lD or quasi-2D exchange 
couplings are highlighted. 

A number of features of these results are of interest. In contrast with the case of F coupled films, bulk and surface 
mode excitations remain well separated in the case of AF coupled layers even in the case where surface exchange 
and anisotropy parameters are identical to bulk values. Differences also arise between these two types of interlayer 
couplings in the case of thin films. In the case of AF coupling, the results are dependent on whether there are an even 
or odd number of layers corresponding to having compensated or un-compensated total moments, respectively. Most 
of our knowledge regarding the sign of exchange interactions is based on rules which may not be entirely applicable 
to surface and thin-film environments.^^ Examination of the spin excitation spectrum in view of these findings offers 
a possible tool to distinguish the type of interlayer coupling. The present work also features a calculation of mode 
amplitudes. These exhibit a variety dependences through the thickness of thin films dependent on the mode in 
question and if the system is quasi-lD or quasi-2D. 

These results generally are applicable to a wide variety of AF coupled layered triangular materials. Among the 
ABX3 compounds, CsMnBrs and CsVCls are notable examples of quasi-lD materials with planar anisotropy.'^ Bulk 
excitations have been well studied in prototypical axial quasi-lD hexagonal materials CsNiCls and CsCoXs.^ Examples 
of quasi-2D compounds with well established bulk magnetic properties include axial CuFeOa^ and planar VX2 and 
LiCrS2.^ Most of the RMnOs materials exhibit planar anisotropy although characterization of thin films in this regard 
is ongoing.- In ultrathin films surface effects dominate and spin reorientation transitions can occur4 A number of 
possible experimental techniques to examine long wavelength surface modes are summarized in our previous workf^i 
Progress using inelastic neutron scattering techniques to probe surface excitation in multilayer systems has been 
reportedfSS Potential extensions of the present calculations might be to include dipolar effects that can be important 
at surfaces even in AF systems. This may require the utilization of purely numerical algorithms such as those based 
on classical equations of motioni^ 
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APPENDIX A 



The definition of matrices A„ and Tn,n' in Eq. (I28|) for a stacked triangular antiferromagnetic thin film with 
easy-axis anisotropy. 



where 
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i3* Cn D 
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y —2 cos(0„ — On') + cr' sin(?„ sin6'„' 



-CT sinW„ sm 
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with matrix elements defined as 

Qn = 2S{l + a)cose,JiO) + 2SD' - S{1 + (7')J'[2-6i^^~6,,^n] 

a„ = S'((l + cr)(cos On + cos^ e„) - sin^ 9n)J{0) - 5(1 - 3 cos^ 0„)D' 

- Sj'{[cOs{6n ~ On~l) +0-'cOs6'„COs6'„_i](l - Sl^n) 

+ [C0s{9n -On+l) +0-' cos 6I„ COS 6'„+i](l -~Sn,N)) 

6„ = -SD'si-a^en 

cf.n = COS 6'„ ± 1 

^tn = COS^ On - {1 + Cr) Sin^ On±l 

= [I- i2S)-']D^. (34) 
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